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The spectral distribution of parametrically excited dipole-exchange magnons in an in-plane mag- 
netized epitaxial film of yttrium-iron garnet was studied by means of frequency- and wavevector- 
resolved Brillouin light scattering spectroscopy. The experiment was performed in a parallel pumping 
geometry where an exciting microwave magnetic field was parallel to the magnetizing field. It was 
found that for both dipolar and exchange spectral areas parallel pumping excites the lowest vol- 
ume magnon modes propagating in the film plane perpendicularly to the magnetization direction. 
In order to interpret the experimental observations, we used a microscopic Heisenberg model that 
includes exchange as well as dipole-dipole interactions to calculate the magnon spectrum and con- 
struct the eigenstates. As proven in our calculations, the observed magnons are characterized by 
having the highest possible ellipticity of precession which suggests the lowest threshold of parametric 
generation. Applying different pumping powers we observe modifications in the magnon spectrum 
that are described theoretically by a softening of the spin stiffness. 



I. INTRODUCTION 

Parallel pumping, the parametric excitation of spin 
waves by means of a linearly polarized microwave mag- 
netic field parallel to an external bias magnetic field, 
was first described in 1960^2 The technique has been 
used for the last five decades to investigate a num- 
ber of interesting phenomena such as the excitation of 
magnctoclastic waves in ferromagncts^ spin-wave para- 
metric instability^^ high density magnon gases and 
condensates^— parametric amplification of spin-wave 
solitons and bulletsj^r— spin pumping and the inverse 
spin Hall effect J2r— While the main principles of para- 
metric pumping have been revealed in experiments with 
bulk ferrimagnetic samples the more recent research is 
focused on the application of this technique to ferro- and 
ferrimagnetic films. Thus, the understanding of the pe- 
culiarities of parametric excitation in the application to 
magnetic films is of paramount importance for the correct 
interpretation of many interesting effects (see, for exam- 
ple, Refs. [H, [HI, and [H focused on microwave driven 
ferrite-platinum bi- layers and nano-contacts) . 

Parallel pumping has two essential advantages in com- 
parison with the direct excitation of spin waves. Firstly, 
when spin waves are directly driven by a microstrip an- 
tenna using an alternating Oersted field, the wavenum- 
ber of the effectively excited spin waves is rather small, in 
the range of 10 3 rad/cm due to the finite width of the an- 
tenna. Therefore, direct excitation of large-wavenumber 
spin waves can hardly be realized. Secondly, parallel 
pumping is a threshold process, i.e. the driving field am- 
plitude h should exceed a certain threshold value ft. C rit, 



which is dependent on the magnetic losses in the mate- 
rial, in order to start the parametric generation. Once 
this threshold has been overcome, the amplitude of the 
parametric spin waves grows exponentially in time un- 
til higher order nonlinear processes limit this amplitude 
to a certain level . 17 ' 18 With this method it is possible 
to achieve much higher spin-wave intensities in compar- 
ison to the direct linear excitation of spin waves with a 
microstrip antenna. 

During the process of parallel parametric amplifica- 
tion, spin waves are excited at half of the initial mi- 
crowave pumping frequency cj p . In terms of energy 
quanta, a microwave photon splits into two magnons 
which are the quasiparticles of the dynamic magnetiza- 
tion. From the laws of the conservation of energy and 
momentum and from the fact that for the case of weakly 
localized pumping the photon wavenumber is negligibly 
small^ it then follows that the pumping process creates 
pairs of magnons with the same frequency lo p /2 and with 
opposite wavevectors k and — k, as shown in Fig. [1] (a). 
Different magnon groups with the same frequency uj p /2 
are parametrically driven at the same time, and only the 
one with the lowest damping survives! 17 ' 20 ' 21 The deter- 
mination of the spectral position of this dominant group 
is of crucial importance not only for the investigation of 
the primarily pumped magnons but for all the magnon 
groups which acquire their energy in subsequent scatter- 
ing processes from the dominant group. 

Here we present a direct observation of the dominant 
group for different pumping frequencies and powers in 
a wide wavevector range up to 1.6 x 10 5 rad/cm using 
Brillouin light scattering (BLS) spectroscopy. The exper- 
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FIG. 1. (Color online) (a) Spectral position of the parametri- 
cally amplified magnon modes using the pumping frequency 
oj p . (b) The elliptical precession of the magnetization M al- 
lows the excitation of magnon modes by an alternating mag- 
netic field h parallel to the external bias magnetic field Hq. 



imcnts arc performed on a single-crystal epitaxial film of 
yttrium- iron garnet (YIG) i 22 i 23 We show experimentally 
that the primarily excited spin waves of the dominant 
group are always located at the 90° branch of the back- 
ward volume magnetostatic wave (BVMSW) having the 
lowest energy. Previous works on this topic were carried 
out by Kabos et al.j&2S. but only for dipolar-dominated 
spin waves with a wavenumber up to 4 x 10 4 rad/cm. 
In our present work we complete the picture by obtain- 
ing more details for the dipolar-dominated waves as well 
as by exploring the parametric processes involving the 
exchange-dominated BVMSW. By changing the pump- 
ing frequency u; p in sufficiently small step sizes, as in- 
dicated in Fig. [T] (a), it is possible to move along the 
dispersion branch of the BVMSW and to determine the 
spectral position of the dominant group. In addition, 
calculations of the magnon dispersion based on an effec- 
tive spin model for YIG have been performed which are 
in excellent agreement with the experimentally obtained 
results. 

This paper is organized as follows: Sec. |TT] provides 
the reader with basic information about the technique of 
parallel parametric amplification. Calculations of the rel- 
evant magnon modes are performed in Sec. Mil In Scc. lIVI 
details of the combined BLS/microwave experiment are 
given, especially the wide wavevector selectivity. Sec. |V] 
shows the experimental results which clarify the position 
of the dominant group for different pumping frequencies 
and powers. These experimental data are supported by 
theoretical investigations of the spin- wave dispersion. Fi- 
nally, in Sec. IVI1 we summarize our results and present 
our conclusions. 



II. PARALLEL PARAMETRIC 
AMPLIFICATION 



the applied field. In this case the component of the mag- 
netization along the external bias magnetic field is con- 
stant: 



Af, 



const. 



(1) 



The parametric amplification can be considered as a re- 
sponse of the spin system to the application of a weak 
alternating magnetic field parallel to the bias magnetic 
field Ho, which is oscillating with a precession frequency 
2wo which is twice the eigenfrequency of the spin wave 
loq we intend to excite parametrically: 



h(t) = ho cos(2cjo*)' 



(2) 



Therefore, for the purely circular precession, the average 
energy pumped into the system vanishes: 



(AE) oc (h ■ M) = (h) M z = h M z 



dt 



cos(2wo0 = 0, 
(3) 

where M denotes the magnetization and To denotes the 
time for one cycle of precession. 

An elliptical precession implies that this circle is 
squeezed in some direction in the plane which is per- 
pendicular to the ^-axis. As shown in Fig. Q] (b), in this 
case the M z component of magnetization oscillates with 
double frequency: 



M z = const. + m Zj o cos(2oj t). 



(4) 



This allows absorption of microwave energy by the spin 
system: 



(AE) oc 



dt 



cos 2 (2u;ot) ^ 0. 



(5) 



For a bulk magnetic material the dynamic demagnetiz- 
ing field along the direction of propagation of the plane 
spin wave gives rise to the ellipticity of the magnetiza- 
tion precession. In this case, the threshold for the parallel 
pumping is given by a simple analytical formula* 1 ^ 



hcrii. 



uj p AH k 
uim sin ?9k 



(0) 



where AH^ = 2Fk/7, Wm = 747rM s , and $k is the angle 
between the propagation direction of the parametrically 
amplified spin waves with wavevector k and the external 
bias magnetic field. Tk denotes the relaxation frequency 
of this magnon mode, 7 is the gyromagnetic ratio, and 
M s is the saturation magnetization. Importantly, from 
this expression it follows that the threshold field h CIlt 
is inversely proportional to the ellipticity £k, which is 
defined by 



A key condition for the parametric amplification is el- 
lipticity of magnetization precession. The circular pre- 
cession implies that the tip of the magnetization vector 
goes around a circle lying in the plane perpendicular to 



£k = 1 



(7) 



Vk 



where m x and m y are the amplitudes of the dynamic 
components of the total magnetization vector M = 



3 



(ma, my, M z ), and y_Lk. Consequently, h CI a reaches the 
smallest possible value at $k = ir/2 when the ellipticity 
reaches its maximum.— ^ 

In the confined geometry of thin films the picture of the 
dynamic demagnetizing field is more complicated. The 
out-of-plane component of the dynamic demagnetizing 
field, caused by the presence of the two film surfaces, 
"compresses" the precession cone in the direction of the 
film plane. As a result, the value of h cr - lt becomes depen- 
dent not only on the direction (as for the bulk materials), 
but also on the value of the in-plane wavenumber. In ad- 
dition, the finite film thickness results in quantization of 
the total magnon wavenumber in the direction perpen- 
dicular to the film surface. A simple analytical expression 
can easily be obtained for the mode of the uniform pre- 
cession (the wave with the vanishing value of the total 
wavenumber) h CI it oc where 



Vk,n = 1 



(8) 



In this work we extend the validity of this formula for the 
entire range of the wavevectors parametrically excited in 
the films. Under this assumption, and with the calcula- 
tion shown in the next section, we find that in the case of 



intermediate film thicknesses (when m Xk 



< 



i) no 



qualitative difference between r] and e exists. Therefore, 
in the following, we restrict ourself to the classical defi- 
nition of ellipticity given by Eq. [7] However, in contrast 
to bulk materials^, in this work the discrete character of 
the magnon spectrum is also taken into account. 

Another particular type of parametric amplification is 
perpendicular pumping. It is characterized by the appli- 
cation of an alternating pumping field in the plane per- 
pendicular to the direction of the equilibrium magneti- 
zation, i.e, the alternating magnetization induced by the 
pumping field is always perpendicular to the external bias 
magnetic field^ In this case, the amplification is caused 
by the parametric decay of the uniform precession mode, 
driven out-of-resonance (i.e. at 2u}q). In our work we use 
the microwave Oersted field of a microstrip transducer 
for the parametric pumping. The Oersted field of such 
a transducer has two components: one in-plane and one 
out-of-plane. The in-plane component is parallel to the 
bias field and thus is responsible for the parallel pump- 
ing process. Similarly, the out-of-plane component of the 
alternating Oersted field is responsible for a possible per- 
pendicular pumping process^ which we want to avoid in 
our work. Fortunately, the in-plane component of the 
Oersted field is maximized below the longitudinal axis 
of the microstrip line and, the out-of-plane component is 
mostly localized near the microstrip edges. Therefore, by 
using a sufficiently small laser spot size and probing the 
area near the longitudinal axis one can exclude the con- 
tribution of the perpendicular pumping mechanism from 
the parametric process. 



III. THEORETICAL INVESTIGATIONS OF 
THE MAGNON SPECTRUM 

In order to find the spectral positions of the paramet- 
rically injected magnons in a film sample, we used a re- 
cently developed microscopic approach^ for calculating 
the magnon energies as well as the ellipticity of the pre- 
cession. Our starting point is a microscopic Hamilto- 
nian that describes the properties of the relevant magnon 
modes in YIG. The model for our magnetic films of a fi- 
nite thickness d contains both exchange and dipole-dipolc 
interactions and is completed by a Zecman term that 
takes into account the external bias magnetic field Hq: 



T~L — — — JijSi ■ Sj — /xHo • Si 



>j 



1 



3(Sj • R,2j)(Sj • R-jj) Sj * Sj 



(9) 



The spin operators Si are normalized such that Sf = 
5(5+1) with an effective total spin 5 per lattice site. The 
sums run over the sites Ri of a cubic lattice with spacing 



12.376 A, and R; 



Hi, /I Riil arc unit vectors in 



the direction of Ry = Ri — Rj. The relevant parame- 
ters are the exchange interactions J,j = J of neighboring 
spins and the magnetic moment [i = g^B, where g is the 
effective g-factor and [Ib is the Bohr magneton. In order 
to proceed using a compact notation we introduce the 
dipolar tensor D°f = D a P(Ri - R,-), 



D 



a/3 



(!-<%)■ 



R^l 3 
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and rewrite our effective Hamiltonian (|9|) as 



I J l0 8 ap 



D 



(10) 



(11) 



ij a/3 



with the z-axis of the frame of reference pointing along 
the direction of the external magnetic field Ho — h/ (ie z . 
We assume that the continuous film is magnetically sat- 
urated by the magnetic field Ho in its plane and the 
ferromagnetic magnetic order is perfect. This allows us 
to expand the Hamiltonian in terms of bosonic opera- 
tors describing fluctuations around the classical ground 
state. Inserting the Holstein-Primakoff transformation^ 
we obtain a bosonized spin Hamiltonian of the form^S 



(12) 



n=2 



Considering the large effective spin 5 = M s a 3 //i ~ 14.2, 
it is sufficient to retain only terms up to n = 2 in the 
formal 1/5 expansion in order to calculate the magnetic 
excitation spectrum. The quadratic part of the Hamilto- 
nian reads 



n 2 



E 



Bi 



bibj 



b]b\ 



(13) 
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with 



i 



B, 



s 



s 



J~\XX 



D vv 



- /)'■ t,d"' i)' 1 ' 1 

2 l^ij ^tj J 



(14a) 



(14b) 



A thin film of thickness d = Na is obviously not transla- 
tional invariant in all spatial directions, such that a full 
Fourier transformation cannot be used. Instead we set 
Hi = (xi,Yi) with r; = (yi,Zi) and use the property of 
the discrete translational invariancc in the y and z di- 
rections to partially diagonalize TL2 via a partial Fourier 
transformation. We expand the bosonic operator 



h = 



1 



^e lk ri 6 k ( It ) , 



(15) 



where N y and N z is the number of lattice sites in the y 
and z direction. The Hamiltonian ()13[) then reads 

H 2 = £]T{[A k ]. y -&]>o& k (*,) 

, [ B k]y 



with the N x N matrices At and B k given by^ 
[A k ]y = 

= SMxij) + S l3 [h + D z 2 



J^e lk r A( Xi -Xj,r) 



(17a) 



[B 



— ik-: 



- Xj,r) 



S 



D?M - 2iD^{ Xij ) - Dl y { Xij )] . (17b) 



The exchange matrix is given by 

- 2(cos(k y a) + cos(fc 2 a))} - S ij+ i - , (18) 

and the dipolar tensor reads 



where £' excludes the term 



yij — Zij 



= when 



0. We can rewrite the quadratic Hamiltonian in matrix 
notation 



Ho 



(20) 



with the grand dynamic matrix 



V k = 



A k B k 



(21) 



and the row vector 

bt=(bl(x 1 ),...,bt(x N )) (22) 
Observing the low symmetry of the matrices 



A k = Ai k A k 
B k = B T k ^ B k 



(23a) 
(23b) 



we have to use a full 2N square transformation 

(24) 



u k -v k 

-wi x k 



to diagonalize the Hamiltonian (|2T))) . This transformation 

connects the vectors of the true magnon operators 7 k = 

(7k l' • • • >7k n) 1° those in the Holstein-Primakoff basis 
via 



U k & k -V k 6l k \ 

-wib k + xW 1 (25) 



k ^ k t ^- k "_ k J 

and transforms to the diagonal Hamiltonian 

% = ^E(7t7„ k ^k(| k k ) + 4 2) (26) 
with the matrix 

fk=(J k t )"X Jy 1 

= Ai&g{E K1 ,...,E KNl E K1 ,...,E KN ) (27) 
and the correction 



4*> = ^TrS k 

k 



(28) 



to the groundstate energy. 



In order to numerically calculate the eigenenergies 
Ek t n, n — 1, . . . , N of all magnon modes and the transfor- 
mation matrix for a given in-plane wavevector k we carry 
out the slowly converging sum in (|19[) using Ewald sum- 
mation technique to set up the grand canonical matrix 
2? k . Following the algorithm of Ref. we first calcu- 
late the Cholesky decomposition /C k of D k , which has the 
property 

= /C k /C k . (29) 
Then we diagonalize the matrix A^ k given by 

In \ 



Mw = /C k / P /C k , I p = 







-1 



JV 



(30) 
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using a unitary transformation U k to obtain the diagonal 
matrix £ k which is connected to the eigenenergies via 



The inverse transformation 



U k W k 

v k x k 



is now obtained from the equation 

r-l 



/C"X El/ 2 . 



(31) 
(32) 

(33) 



Having constructed the full transformation, we then ob- 
tain the amplitude structure of the modes across the film. 
It is worth noting that in our approach the anti-nodal 
points of the magnon modes are located on the film sur- 
faces which is consistent with the unpinned surface spin 
configuration.— 

In line with the discussion in the previous section we 
also calculate the ellipticity of the modes e ki „. The el- 
lipticity can be calculated from the expectation values of 
the spin operators in the magnon eigenstates^ 



1 



<(^) 2 >k,„' 



(34) 



S a denotes the a component of the total spin operator 
and the expectation value has to be taken in the magnon 
cigenstate |k, n) = 7 kn |0) of the n-th magnon mode. 
According to a calculation in lowest order 1/5 we use 
the expansion of the Holstein-Primakoff transformation 
to express the spin operators in terms of boson operators 
to obtain 



({s x )\ 



s 



[(fr-k ■ &k)k,« + (^_ k ' &_ k )k,n 
+ (bl ■ $- k W,n + (bt ■ 6 k )k,n] 



<(5 y ) 2 ) k ,„ - - [<L k • & k )k,n - (L k • &L k )k,n 



+ (bl-bl k )^ 



(35a) 



(35b) 



An inversion of the transformation (|25[) allows the cal- 
culation of the expectation values, 



(6- 



(bt ■ 6 k ) k , ra 



3 
3 

(A ■ 5L k >M = YKl U t}nAV k } jn + [Wl\ nj [X k ] jn ), 
3 

<&-k • bk)k,„ = ^(WXiW^Jn + [xtU[w k } jn 



(36a) 
(36b) 



(36c) 



(36d) 
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FIG. 2. Theoretical result for the spectrum with M B = 1742 G 
(a) and the ellipticity of the corresponding magnon modes 
as a function of the wavevector (b) for a YIG film of thick- 
ness d = 5 j-im in a magnetic field of Ho = 1750 Oe. 
Our theoretical approach delivers the spectral position of all 
magnon modes for the model including both dipole-dipole 
and exchange interactions. The steepest section of the dis- 
persion curves at small wavenumbers are the surface Damon- 
Eschbach mode which is strongly hybridized with the higher- 
order thickness modes (see e.g. Refs. l34l and l35h. Note that 
the ellipticity ek,n decreases as a function of the wavevector 
and with growing mode number n. 



such that the ellipticity is given by 

2£>k.n 



fk.. 



(37) 



, n 

with the abbreviations 

3 

+ \[U k ] jn \ 2 + \[W k ] jn \ 2 ), (38a) 
B Kn = 2 Re XX[Vk] jn [U£\ Jn + [X k ] jn \W£] jn ) . (38b) 



In the following we only analyze the characteristics 
of magnons with d k = 7r/2 as these magnons have the 
highest ellipticity and consequently one may expect the 
strongest parametric coupling to this type of magnons. 
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Furthermore, in our experiment these magnons propa- 
gate along the length of the microstrip transducer. This 
length is much greater than the free propagation path for 
the magnons. Therefore, one may expect that, in con- 
trast to all other magnons which have a non-vanishing 
component of the in-plane wavevector perpendicular to 
the longitudinal axis of the transducer, the threshold of 
the parallel pumping excitation for them is not affected 
by the energy loss due to magnons escaping the pumped 
area. 27 ' 36 Figure [5] shows the calculated magnon disper- 
sion along with the ellipticity of the precession as a func- 
tion of the magnon wavenumber k. The spectrum com- 
prises the finite set (limited by the number of crystallo- 
graphic elementary cells of 12.376 A on the film thick- 
ness) of the quantized volume magnon modes (BVMSW) 
and a steep branch corresponding to the surface Damon- 
Eshbach magnon mode. 



IV. SETUP 

The measurements were performed using combined mi- 
crowave and optical facilities where the magnon system 
was parametrically pumped by the microwave circuit, 
and the response of the magnon system is analyzed by 
means of Brillouin light scattering spectroscopy. The mi- 
crowave circuit comprises a microwave source, a switch, 
and a microwave amplifier connected to a probe section 
(see Fig. [3] (a)). In contrast to the conventional approach 
where the pumping frequency is held constant and the 
applied bias magnetic field Hq is swept we examined the 
magnon spectrum at different pumping frequencies hold- 
ing the applied field and consequently the ground state 
for small amplitude-excitations constant. This implies 
that the spectral characteristics shown in Fig. 2 were pre- 
served during the experiment. The pumping microwave 
Oersted field was created by a 50 /im wide short-circuited 
microstrip line. The regular microstrip line was utilized 
instead of the commonly used microwave resonato r 8 i 13 i 27 
as it enables us to change the pumping frequency without 
the complicated realignment of the probe section. The 
investigated sample, a 15 mm long and 3 mm wide YIG 
film with a thickness of 5 ^m, was placed on top of the 
microstrip line. An external static field Hq of 1750 Oe 
was in the film plane and perpendicular to the longitu- 
dinal axis of the microstrip. In order to excite dipolar 
as well as exchange-dominated magnons, the frequency 
of the microwave source was varied from 13.6 GHz up 
to 14.6 GHz in steps of 20 MHz, which corresponds to 
change lu p /2 from 6.8 GHz to 7.3 GHz in steps of 10 MHz. 
Thus, according to the magnon spectrum, a wide range 
of wavenumbcrs from zero up to 1.6 • 10 5 rad/cm could 
be investigated by means of BLS spectroscopy. 

The experiments were performed at different mi- 
crowave powers from 100 mW to 10 W. In order to reduce 
any possible thermal effects which might have poten- 
tially influenced the magnon spectrum at high pumping 
powers by modifying the saturation magnetization M s , 




FIG. 3. (Color online) (a) Geometry and block diagram for 
the experimental setup of the wave- vector selective BLS. (b) 
Geometry for the scattering of the probing light on the Bragg- 
grating generated by spin waves in the BLS technique. 



anisotropy fields, and the exchange stiffness constant the 
pumping was applied in 2 /is microwave pulses separated 
by 20 fj,s time intervals. For the same reason, a metal- 
lized aluminum nitride substrate, which is known for its 
high thermal conductivity, was chosen as a base plate for 
the probe section. 

In order to probe the parametrically pumped magnons, 
the BLS measurements were performed in the back scat- 
tering geometry, where a single-mode solid state laser 
with a wavelength of 532 nm was focused onto the sam- 
ple by using an objective with a high numerical aper- 
ture. The focal point of 20 /jm in diameter was centered 
on the longitudinal axis of the microstrip line where the 
microwave Oersted field is parallel to the bias magnetic 
field Ho, and thus the condition for parallel parametric 
pumping is satisfied. The backscattered light was col- 
lected with the same objective and sent to a tandem 
Fabry-Perot interferometer for frequency and intensity 
analysis. In a classical description, Brillouin light scat- 
tering can be interpreted as the diffraction of the probing 
light from a moving Bragg grating produced by a magnon 
mode, see Fig. [3] (b). As a result, some portion of the 
scattered light is shifted in frequency by the frequency 
of this mode (in our case ±(J p /2). The intensity of the 
inelastically scattered light is directly proportional to the 
number of magnons in the mode. In addition, the diffrac- 
tion from the grating leads to a transfer of momentum 
during this process. By changing the angle Ob between 
the sample and the incident light beam, the wavevector 
selection with a resolution of up to 4500 rad/cm in a 
film plane can be implemented as shown in Ref. l37l The 
incident angle Ob of the probing light determines the 
selected magnon wavenumber k = 2fc sc sin(©B), where 
the wavenumber of the scattered light k sc is equal to the 
wavenumber of the incident light k- ln . 



V. EXPERIMENTAL RESULTS 

Figure @] shows the measured BLS intensity of para- 
metrically injected magnons as a function of the fre- 
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quency and the wavevector for different pumping pow- 
ers. To increase the signal-to-noise ratio the intensity 
was integrated across the entire width of the inelasti- 
cally scattered peak for w p /2. The color-coded intensity 
maps, where blue (dark) corresponds to low intensities 
and orange (bright) to higher intensities, were recorded 
by changing the pumping frequency for the given inci- 
dent angle 0b of the probing light (see Fig. [1} . Once the 
frequency measurement cycle had been completed, the in- 
cident angle Ob w &s changed and the measurement were 
repeated in this way for wavenumbers ranging from to 
1.6 x 10 5 rad/cm in steps of 8.2 x 10 3 rad/cm. 

In Fig. [4] (a), the color coded BLS-intensity map mea- 
sured at the pumping power Pi =0.1 W, which is slightly 
(approximately 1 dB) above the threshold power of the 
parametric generation, is compared with the theoretical 
calculations of the relevant modes of the magnon spec- 
trum. It is clearly visible that the detected magnons are 
located along the dispersion curve representing the low- 
est magnon mode n = 1 while other regions of the spec- 
trum show practically no increase in the BLS intensity. 
This behavior can be attributed to the fact that from all 
magnon groups pumped at the same time only the dom- 
inant group having the lowest damping and the highest 
coupling to the pumping field is significantly populated. 

The lowest mode corresponds to the backward volume 
magnetostatic wave having one node along the thickness 
of the YIG film. Modes with higher frequencies belong 
to higher order standing-wave modes, quantized perpen- 
dicular to the film plane. As we have no reason to as- 
sume different damping for different modes spread across 
the relatively narrow frequency and wavenumber ranges, 
our attention must be focused on the coupling efficiency. 
Since the elliptical precession of the magnetization allows 
one to excite the corresponding magnon modes (compare 
Fig. [TJ, as discussed above, we assume that the thresh- 
old of parametric generation is proportional to the in- 
verse ellipticity, ft, cr ;t oc l/ek,n- Observing Fig. [SJ one 
notices that the lowest magnon mode, n = 1, has the 
largest ellipticity, and correspondingly the lowest thresh- 
old, in the entire pumping frequency range. Therefore, 
one would expect parametric amplification of this mode 
as the dominant magnon mode. 

While the bias magnetic field is kept fixed during the 
measurement, the saturation magnetization M s is in- 
fluenced by the pumping of magnons. Further, single 
magnons carry a total spin of 1 and every magnon re- 
duces the magnetization by one Bohr magneton In 
order to account for that in the theoretical calculation 
of the spectrum we fix the exchange coupling J at zero 
pumping and set the relevant spin stiffness to 



Since the absolute magnon densities, composed by ther- 
mal magnons and pumped magnons, are unknown at this 
point, we use the saturation magnetization as a fit pa- 
rameter. An adjustment of M s then changes the posi- 




Wavenumber (*10 4 rad/cm) 

FIG. 4. (Color online) Color-coded BLS-intensity map for the 
three different pumping powers of 10 W, 1 W and 100 mW. 
In panel (a) the white solid curves show the theoretically cal- 
culated dipolar-exchange spectrum of transversal spin waves 
propagating across the magnetization direction. In panels (b) 
and (c) the single white solid curves correspond to calculations 
of the lowest transversal magnon mode. The saturation mag- 
netization has been used as a fit parameter: (a) M s = 1742 G, 
(b) M s = 1715 G, (c) M s = 1676 G. 



tion of the ferromagnetic resonance and the slope of the 
spectrum that originates from the effective exchange cou- 
plings. However, the slope of the spectrum also depends 
on the saturation magnetization and therefore needs to 
be adjusted to the experimental data. The solid curves 
in Fig. IUJb)-(c) correspond to the transversal BVMSW 
mode having the lowest energy. To fit the experimental 
data we varied the value of the saturation magnetiza- 
tion. We saw that with increasing pumping power, M s 
decreased significantly, in accordance with the theory, 
from 1750 G for the undisturbed sample to 1676 G at 
pumping level of P3 = 10 W. At the same time, one can 
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FIG. 5. (Color online) Inverse ellipticity curves calculated as 
a function of frequency for the magnon modes propagating 
across the magnetization direction in a YIG film of thickness 
d — 5 /im. The saturation magnetization M s = 1742 G cor- 
responds to the pumping power of Pi = 0.1 W. Since the 
ellipticity ek.n is largest and the corresponding threshold of 
parametric generation /i cr it oc l/ek,n is lowest for the lowest 
magnon mode, n = 1, pumping to this mode is expected at 
any given frequency. 

see that the precision of the fit decreases for the largest 
pumping power. It can be connected to a spread of 
magnon population beyond the dominant magnon mode 
and to corresponding blurring of the BLS intensity over 
the wide frequency and £;-number ranges. This spread is 
caused by a number of physical phenomena. Firstly, at 
high pumping levels a number of frequency-degenerated 
magnon groups with damping higher than one for the 
dominant group can be excited^ Secondly, two-magnon 
scattering processes from crystal defects, dopants, and 
other static magnetic non-uniformities breaks the law of 
conservation of momentum, and thus lead to a spreading 
of the parametrically excited magnon population along 
the isofrcquency lines in the k-spacc. Thirdly, interac- 
tions between magnons and lattice vibrations will cause 
a redistribution of the magnons towards states with lower 
energy.— Finally, non-elastic four-magnon scattering pro- 
cesses are responsible for occupation of energy levels 
around the initial magnon state at half of the pumping 
frequency, and thus for consequent thermalization of the 
parametrically injected magnons across the entire spin- 
wave spectrum.— 

In addition, we want also to comment on the role of 
the intermodal concurrency in parametric generation^ 
We have observed no fine structure in the occupation 
of quantized standing modes at the beginning of the 
magnon spectrum as it was predicted by the theory in 
Refs. HH and and reported in Ref. 0- This might be 
due to the fact that under our experimental conditions 



the inverse ellipticity curve of the lowest BVMSW mode 
lies below the corresponding curves of the highest stand- 
ing modes and, as a result, this mode is solely amplified 
across the entire range of the experimentally probed fre- 
quencies (Fig. [5]). Moreover, in accordance with our es- 
timation based on Eq. [37l the inverse ellipticity curve of 
the lowest BVMSW mode intersects first with the 65th 
standing mode at 9.6 GHz. The in-plane wavenumber of 
the lowest mode at this frequency is k = 4 x 10 5 rad/cm. 
The wavenumber of the 65th mode, calculated across the 
film, is nir/d = 4.08 x 10 5 rad/cm. Both of these values 
are significantly larger than the highest optically acces- 
sible wavenumber value of 2.6 x 10 5 rad/cm. As a result 
no intermodal jumps in the magnon generation, and con- 
sequently no discrete structure in the BLS intensity map 
are expected to be observed in our experiment. 

VI. CONCLUSION 

We have observed the generation of parametric 
magnons under parallel pumping using the wavevector 
resolved BLS technique. In the theoretical part we con- 
structed the magnon eigenstatcs for a thin ferromag- 
netic film in the framework of a microscopic Heisenberg 
model in linear spin-wave theory to be able to calculate 
such physical characteristics as eigenfrequencies, elliptic- 
ity, and spatial distribution of magnon modes. A good 
agreement between the experimentally determined spec- 
tral position of photon-coupled magnon pairs in a tan- 
gcntially magnetized YIG film and the lowest-frequency 
magnon mode propagating in a film plane perpendicu- 
larly to the magnetization direction was obtained. Com- 
bining the theoretical results with the experimental ob- 
servations, the dominant parametric excitation of the 
lowest mode is understood as a result of its highest el- 
lipticity in the range of accessible wavenumbcrs. A fine 
structure which was previously reported in the threshold 
of parametric generation in YIG films was not detected in 
our experiment. This is possibly caused by the structure 
of the magnon spectra of the magnetic films in our case: 
the ellipticity of the lowest-frequency transversal magnon 
mode may be maximal across the entire range of the op- 
tically detectable fc-numbers. In addition, a significant 
broadening of the populated spectral area is observed 
at high pumping powers. This phenomenon is associ- 
ated with a variety of non-linear and linear scattering 
processes leading to thermalization of the parametrically 
injected magnons. 
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